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Stars with masses in the range 7—10 M & end their lives either as massive 
white dwarfs or weak type II supernovae, and there are only limited obser- 
^ ! vational constraints of either channel. Here we report the detection of two 

O white dwarfs with large photospheric oxygen abundances, implying that they 

are bare oxygen-neon cores and that they may have descended from the most 
massive progenitors that avoid core-collapse. 

White dwarfs represent the endpoints of stellar evolution for the overwhelming majority of 
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all stars. Most white dwarfs in the Galaxy have carbon-oxygen (CO) core compositions, be- 
ing descendants of low to intermediate mass stars that underwent hydrogen and helium core 

ON ' 

burning. Stars with initial masses 7M & < M < 1OM will reach sufficiently high core temper- 
atures to proceed to carbon burning, and produce either oxygen-neon (ONe) core white dwarfs, 
or undergo a core-collapse supernova (SN II) via electron capture on the products of carbon 
burning (0). The exact outcome of stellar evolution in this mass range depends critically on the 
detailed understanding of mass loss ©, the relevant nuclear reaction rates and the efficiency 
of convective mixing in the stellar cores ©. Some observational constraints on stellar mod- 
els come from analysis of SNII progenitors ©, which suggest a lower limit on the progenitor 
masses of 8~t\'®M Q . A linear extrapolation of the observed relation between the masses of white 
dwarfs and their progenitors © up to the maximum mass of white dwarfs leads to a broadly 
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consistent result. 

Information about the core compositions of white dwarfs, in particular of those descending 
from the highest possible masses, would help to test and improve the theory of stellar evolution. 
Unfortunately, almost all white dwarfs have hydrogen and/or helium envelopes that, while low 
in mass, are sufficiently thick to shield the core from direct view. Asteroseismology has the 
potential to unveil their core composition ©, but observational studies attempting to exploit 
this potential remain ambiguous [see e.g. the case of GD358, (0IHI21)]. Core material can be 
directly detected in the photospheres of a small number of stars (9) that underwent a very late 
thermal pulse during their asymptotic giant branch (AGB) evolution, ejecting a large fraction 
of the envelope (|7Q|) and leaving a white dwarf with only a thin layer of helium. Examples of 
hydrogen-deficient post-AGB stars are the hot PG1159 objects, white dwarfs with helium-rich 
(DB) atmospheres and cool carbon-rich (DQ) white dwarfs. The most extreme case known 
to date is 1H1504+65, which is hydrogen and helium deficient, representing a hot, naked CO 
surface (fTT]) . 

Recently, the Sloan Digital Sky Survey (SDSS), through its comprehensive spectroscopic 
snapshot of the Galactic stellar population, has revealed a small class of hydrogen-deficient 
white dwarfs (|72l) . Their spectra are consistent with nearly pure carbon atmospheres (|73)) . It 
has been suggested that these "hot" DQ white dwarfs represent the evolutionary link between 
objects such as PG1 159 and 1H1504+65, the DB white dwarfs, and the cool DQ white dwarfs 
(23123). 

Until now, all known white dwarfs for which photospheric oxygen and carbon abundances 
have been determined have abundance ratios O/C < 1. Some stellar models predict that the 
most massive stars avoiding core-collapse will result in ONe white dwarfs with very low carbon 
abundances (Ti6l) . Should such a core lose its hydrogen envelope, an extremely oxygen-rich 
spectrum would be expected. Here we present the results of a search for white dwarfs with large 
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photospheric oxygen abundances within the SDSS spectroscopic Data Release (DR) 7 (|77b . 

We selected all spectroscopic objects in DR7 that fall within the ugriz color space of the 
white dwarfs from (|78l) . We then subjected these spectra to an automatic measurement of the 
equivalent widths [see ([791) for details of the method] of the O I 616 nm and 778 nm multiplets, 
which are the strongest O I lines expected in a cool, oxygen-rich atmosphere. Another strong 
O I multiplet is located at 845 nm, however, this region can be affected by residuals from the 
night sky line subtraction, and we therefore did not include it in our procedure. As an additional 
constraint, we required the spectra to have a signal-to-noise ratio > 10 in the regions around the 
two O I lines, resulting in the analysis of 25639 SDSS spectra (|20l) . 

We found 1000 spectra with a formal > 4cr detection of both 01 lines, and visual in- 
spection of these spectra showed that 998 were erroneously flagged by our automated rou- 
tine because of poor data quality. Only two objects with genuine OI 616 nm and 778 nm ab- 
sorption, SDSS J092208. 19+2928 10.9 and SDSS Jl 10239.69+205439.4, survived our scrutiny. 
Both have high proper motions, 0.275'7yr and 0.163'7yr, respectively (1271) . consistent with them 
being nearby low-luminosity blue stars. Aside from the 1 616 nm and 778 nm absorption lines, 
the spectrum of SDSS 0922+2928 (Fig. 1) is typical of a cool DQ white dwarf, with strong 
C 2 Swan bands in the blue, superimposed by a number of atomic carbon lines. In contrast 
SDSS 1 102+2054 is a unique white dwarf (122)) with a photospheric spectrum totally dominated 
by absorption lines of OI (Fig. 1). The only other noticeable features are some weak lines of 
atomic C I as well as Ha and H/5. 

We determined the effective temperatures of the two white dwarfs by fitting their SDSS 
ugriz photometry with two grids of theoretical atmosphere models: one for a pure He com- 
position [Supporting Online Material (SOM), Sect. 1], and one including H, C, and O as well 
(SOM, Sect. 2). The white dwarf temperatures and chemical abundances determined from these 
fits depend only mildly on the assumed surface gravity log g and, because their masses are not 
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constrained by the available observations, we fixed M = O.6M , corresponding to the canoni- 
cal value of log g = 8 [ @2J, SOM Sect. 2]. 

Independent of the details of the atmospheric composition, the photometry suggests that 
both white dwarfs have temperatures in the range 8000-9000 K (Table 2). The implied distances 
depend strongly on the assumed surface gravity via the radius, here we used the mass-radius 
relation of (|2?1) . We used the effective temperatures (and assumed surface gravity) to calculate 
synthetic spectra for a final determination of atmospheric abundances from the absorption lines 
in the SDSS spectroscopy (SOM, Sect. 2). 

For SDSS 1 102+2054, we adopted T cff = 8150 K and log g = 8 and iterated the abundances 
of H, C, and O, until we obtained a reasonable fit to the line strengths. Because this was not 
possible for all oxygen lines we concentrated on the strong lines at 778 and 845 nm, which 
originate from the lowest lying levels of all optical lines. These lines depend only weakly on 
temperature and are best suited for a determination of abundances. The best compromise was 
achieved for the abundances in Table 3. 

However, for this model the carbon multiplet near 712 nm and the 834 nm line, Ha, and 
many OI lines from higher levels are much too broad, indicating too high a pressure in the 
atmosphere. In addition, O I lines from the higher levels are much too weak. The most likely 
solution of this discrepancy is to assume a higher temperature, which increases the strength 
of the high-excitation lines compared to the strong multiplets (SOM, Sect. 2). A reasonable 
compromise was achieved with T cff = 10500 K, log g — 8 (Table 3). The remaining differences 
in the absorption lines between the model and observation are likely within the uncertainties of 
the atomic data and in particular the van der Waals broadening constants (SOM, Sect. 2). 

This model, however, is in significant disagreement with the temperature derived from the 
SDSS photometry. One plausible cause for the discrepancy between the temperatures deter- 
mined from our fits to the photometry and to the spectroscopy are uncertainties in the ultraviolet 
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opacities of an atmosphere of such peculiar composition, which strongly affect the slope of the 
optical continuum [ d2?]), SOM Sect. 1]. 

The atmospheric parameters of SDSS 1102+2054 thus remain somewhat uncertain. Never- 
theless, when varying the parameters within the range demanded by photometry, line strengths 
and widths, the C and O abundances did not change more than factors of ~ 3. Thus the O/C 
abundance ratio in SDSS 1 102+2054 is clearly much larger than one, which makes this star 
unique among the many thousands of known white dwarfs. 

Because, according to the photometry, the effective temperature for the two white dwarfs 
cannot be very different, the carbon abundance in SDSS 0922+2928 must be higher in order 
to produce the strong C2 Swan bands. On the other hand, the hydrogen abundance must be 
lower, because no Balmer lines are visible. We fixed the H/He ratio at 10~ 5 , the upper limit 
allowed by H/5. Further reduction did not substantially influence the models because the major 
electron donors are C and O. Using T c g = 8270 K from the fit to the photometry, and assuming 
logg = 8.0, we obtained a good fit to the spectrum with the oxygen abundance decreased by 
~ 0.2 dex and a carbon abundance increased by ~ 0.6 dex with respect to SDSS 1102+2054 
(Table 3). Therefore, also in SDSS 0922+2928, O/C > 1. With the exception of the strong O I 
multiplets 778 and 845 nm most of the higher excitation O I lines and the C I lines show the same 
problem as in SDSS 1 102+2054: they are too weak and too broad. This is again much improved 
at a hotter temperature, e.g. 10000 K, which is, however, in conflict with the photometry (SOM, 
Sect. 1). 

The low H abundance suggests that SDSS 0922+2928 and SDSS 1 102+2054, similar as the 
other classes of hydrogen-deficient white dwarfs, underwent a late shell-flash, leaving a He- 
dominated atmosphere. At the low temperatures of these stars, the He convection zone extends 
sufficiently deep to dredge core material up into the atmospheres (|25II2(5I) . Given the large age of 
the two stars (> 500 Myr), gravitational diffusion will unavoidably lead to a larger concentration 
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of carbon in the envelope. The only plausible explanation for the observed O/C abundances is 
that these two white dwarfs have overall very low carbon mass fractions, and hence represent 
naked ONe cores. As such, they are distinct from 1H1504+65 and the "hot" DQ white dwarfs. 

Most stellar models that produce ONe cores predict a layer of CO surrounding the core that 
should be sufficiently thick to avoid upward diffusion of large amounts of oxygen. However, 
a sequence of stellar evolution calculations approaching the mass boundary of stars forming 
ONe white dwarfs and those undergoing electron-capture SN II show that the mass of the CO 
layer decreases with increasing initial stellar mass (|27][2#][76]). Thus SDSS 0922+2928 and 
SDSS 1102+2054 may have descended from the most massive stars avoiding core-collapse, in 
which case they would be expected to be very massive themselves. Our current data is insuffi- 
cient to provide any unambiguous measure of the masses of SDSS 0922+2928 and SDSS 1 102+2054, 
however their C/He abundance is at the top of the range observed in the carbon-rich sequence 
of "cool" DQ (|23ll29l white dwarfs, some of which have parallax measurements implying that 
they are of high mass. 

We explored spectral models for SDSS 1102+2054 with the abundances kept fixed as in 
Table 3, but adopting surface gravities of logg = 8.5 and 9.0, which correspond to masses of 
0.9 M Q and 1.2M Q (SOM, Sect. 2). Broadly similar fits to the absorption lines can be achieved 
for higher surface gravities if the temperature is increased by 1000 to 2000 K as well. For 
log g = 9.0, the strongest O I lines become somewhat too broad compared to the observations, 
and we conclude that the currently available data is consistent with masses of up to ~ 1M . 

Initial models of the evolution of intermediate mass stars predict that ONe cores should 
also contain significant amounts of magnesium. Updated nuclear reaction rates have lead to a 
substantial downward revision of the Mg abundances (T3Q|) . For SDSS 1 102+2054, which has 
the better quality spectrum of the two white dwarfs presented here, we can place an upper limit 
on the magnesium abundance of log[Mg/He] < —6.1 from the absence of the Mgll 448 nm 
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line. 
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Fig. 1. Sloan Digital Sky Survey spectroscopy. The observed spectra of SDSS 0922+2928 
(bottom) and SDSS 1 102+2054 (top) are shown as thick gray lines. Superimposed as thin black 
lines are models for T cff = 10000 K and log 5 = 8 (bottom) and T cff = 10500 K and logg = 8 
(top). The wavelengths of the strongest OI and CI absorption lines, as well as of Ha and H/3 
are indicated by tick marks. 
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Table 1. Equatorial coordinates, magnitudes and proper motions. 





SDSS 0922+2928 


SDSS 1102+2054 


RA(2000) 


09 22 08.19 


11 02 39.69 


Dec(2000) 


+29 28 10.9 


+20 54 39.4 


u 


18.77 ±0.02 


17.58 ±0.02 


9 


18.52 ±0.02 


17.24 ±0.01 


r 


18.47 ±0.02 


17.26 ±0.01 


i 


18.58 ±0.02 


17.34 ±0.01 


z 


18.71 ±0.03 


17.46 ±0.02 


p.m. 


0.275" yr- 1 


0.163" yr" 1 



Table 2. Atmospheric parameters for a surface gravity of log g = 8. 



System 


T e s 


[K] 


distance [pc] 




pure He 


with metals 


pure He 


with metals 


SDSS 0922+2928 


8720 ± 260 


8270 ± 320 


141 ±6 


122 ±7 


SDSS 1102+2054 


8820 ± 110 


8150 ± 150 


81 ±2 


67 ±2 



Table 3. Photospheric abundances (by number). 



System (T eS /logg) 


log(H/He) 


log(C/He) 


log(0/He) 


log(0/C) 


SDSS 1102+2054 (8150/8.0) 


-3.2 


-3.2 


-1.8 


1.4 


SDSS 1102+2054 (10500/8.0) 


-4.0 


-3.6 


-1.8 


1.8 


SDSS 0922+2928 (8270/8.0) 


-5.0 


-2.6 


-2.0 


0.6 
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Supporting Online Material for 
Two white dwarfs with oxygen-rich atmospheres 
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We describe here in more detail the white dwarf model atmospheres that we have computed 
for the analysis of the SDSS data of SDSS 0922+2928 and SDSS 1 102+2054, and outline the 
approach of our spectral fits. Two grids of atmosphere models were calculated. The first grid 
consisted of pure helium models, used to provide a first estimate of the effective temperatures 
of both white dwarfs from fitting their SDSS photometry. The second grid, including H, C, and 
O, was used to refine the temperature estimate and to determine the atmospheric abundances 
from fitting the observed absorption lines in the SDSS spectra. 

1 He-model fits to the SDSS photometry 

In order to determine a first estimate of the effective temperatures of the two white dwarfs, we 
fitted the He-models to the observed ugriz magnitudes, which were corrected to the AB scale by 
using the corrections suggested by the SDSS project (-0.04, 0.00, 0.00, 0.01, 0.02^. Because the 
mass of the two white dwarfs is currently unconstrained by observations, we fixed the surface 
gravity to the canonical value of log g — 8.0. Good fits were achieved for both objects, implying 
temperatures in the range 8000-9000 K (Table 2, Fig. SI). We note that for both objects the 
SDSS data exhibits a discrepancy between the u-band flux implied by the photometry, and the 
extrapolation of the flux-calibrated spectroscopy. As part of our search for white dwarfs with 
oxygen-rich white dwarfs, we have inspected the SDSS data of a large number of objects, and 
find that such a discrepancy is rare, but not truly exceptional, and based on the available data, 
it is not possible to give preference to either the photometry or spectroscopy. Taking the two 
spectral energy disributions at face value, a fit to the photometry is bound to result in a lower 
temperature compared to a fit to the spectroscopy. We will come back to this point below. 

1 http://www.sdss. org/dr7/algorithms/fluxcal.html#sdss2ab 
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2 Fitting the SDSS spectroscopy with variable H,C, and O 
abundances 



Since He is practically neutral at the temperatures estimated above, the three other elements 
(H, C, O) provide a significant electron density, which correspondingly changes the opacity 
and thermal structure of the atmospheres. The equation of state included 6 diatomic molecules 
(H 2 , CH, OH, C2, CO, O2). In particular the very stable molecule CO is important, because it 
reduces the atomic carbon abundance considerably. Also included in the second grid was the 
continuous absorption of C and O using cross sections obtained from the TOPBASE database 
of the Opacity Project ([DE])- The necessary O I and C I line data were obtained from the VALD 
database (0) and from the Kurucz list (CDROM23). Many of the 01 lines from Kurucz list 
that were prominent in the model were not visible or much weaker in the observed spectra, 
indicating that some of this line data is wrong. We therefore gave preference to the VALD data, 
where the problematic lines from the Kurucz list are either absent, or have gf values several 
orders of magnitude smaller. Recent calculations (4) confirm the data for atomic oxygen in the 
databases we used. 

There are a number of 01 lines where the core electrons are not in the ground state A S° 
configuration, but in the 2 D° or 2 P° state. In these cases the core electrons (without the excited 
electron which makes the transition) are not in the configuration of the O II ground state and 
the ionization therefore would occur to an excited level of the O II ion. This has to be taken 
into account in our approximate calculation of the van der Waals broadening constants, which 
depend on the distance of the levels to the ionization limit. We use the ionization limit of that 
series, that is 16.942 or 18.636 eV instead of the normal 13.618 eV. From that limit and the 
term levels we estimate the broadening constant using standard approximations (5). Similar 
approximations are also used for the Stark broadening constants. 

We re-fitted the SDSS ugriz photometry with H, C, and O abundances close to the final 
values, which resulted in slightly lower temperatures for the two stars compared to the pure 
He-models (Table 2). In a final step, we varied the abundances of H, C, and O while keeping the 
effective temperatures fixed (8270 K for SDSS 0928+2928 and 8150 K for SDSS 1102+2054) 
and assuming log g = 8.0. In this process, we focussed on the strong OI lines at 778 and 
845 nm, which originate from the lowest lying levels of all optical lines and which only mildly 
depend on the prevailing effective temperature. It turned out that no abundance pattern would 
provide a fully satisfying fit to all observed lines, and we consider the abundances listed in 
Table 3 to represent the best compromise. 

A specific problem with this model is that the O I lines from the higher levels are much too 
weak. One way to reduce the pressure in the atmosphere models is to reduce the amount of He. 
We have calculated a series of models between 7000 and 1 1000 K with the O/He ratio increased 
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up to values > 10. In all cases the 01 lines are much too strong, and we rule out this possi- 
bility. The most likely solution to this problem is then to assume a higher temperature, which 
increases the strength of the high-excitation lines compared to the strong multiplets. Assuming, 
as before, a canonical log g — 8, a reasonably good fit to the absorption lines was achieved 
with T e fj = 10500 K (Fig. 1). However, the slope of this model is substantially bluer than the 
SDSS ugriz photometry (Fig. SI). A fundamental uncertainty in modelling white dwarf spec- 
tra with such unusual abundances is our limited knowledge on the sources of ultraviolet opacity 
in their atmospheres. Any flux absorbed at ultraviolet wavelengths will be re-distributed to 
longer wavelengths, and hence result in a steepening of the slope in the optical bandpass. An 
example of the marked differences in the amount of ultraviolet opacity between pure He and C 
atmospheres has been illustrated by ©. Our models predict strong ultraviolet absorption edges 
of O I and C I (which we include with the values obtained from TOPBASE), and a strong Let 
line broadened mostly by neutral helium (which we describe with the theory developed by our 
group). The currently available data does not provide quantative constraints on the ultraviolet 
characteristics of SDSS 0922+2928 and SDSS 1102+2054, and hence a systematic uncertainty 
on their effective temperatures remains. We conclude that, within the framework of the spec- 
troscopic analysis, it is reasonable to prefer the spectroscopic solution of the higher effective 
temperatures, as it is impossible to satisfactorily model the observed ratio of O I lines with high 
and intermediate excitation levels at the low (photometric) temperature. 

As discussed in the main manuscript, the large O/C abundance ratios found for these two 
white dwarfs strongly suggests that they represent bare oxygen-neon cores that may have de- 
scended from massive progenitor stars. Under this assumption, their masses based on any of 
the current initial mass-final mass relations ([71 I^IPl liOl 1771 1721) are expected to be ~ 1 M , cor- 
responding to a likely range for the surface gravity of 8.5 < log g < 9.0. In order to explore the 
effect of a higher gravity on the spectral fit, we have calculated a grid of models spanning the 
range log g = 8.0 to 9.0. Since the spectrum is largely dominated by the ionization balance of 
He and the other elements, a change of log g can to some extent be compensated by a change 
in temperature, and to first order the abundances do not change significantly. In fact, we do 
find that broadly similar fits can be achieved when increasing T c ^ along with \ogg, but keeping 
the abundances fixed at values given in Table 3. Figure S2 illustrates the fits with our canonical 
model (T c fj = 10 500 K, log g = 8.0), and models for higher surface gravities, (T e fj = 11 500 K, 
log g = 8.5) and (T cff = 12 500 K, log g = 9.0). For log g — 9.0, the strongest O I lines become 
somewhat too broad compared to the observations. We conclude that the currently available data 
are consistent with a broad range of white dwarf masses, including relatively massive (~ 1 M Q ) 
stars, and that the main conclusion of a large O/C abundance ratio is independent of the mass 
assumed for the spectral fitting. 
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Supporting Figure 1. SDSS spectroscopy (gray lines) and photometry (green dots) of SDSS 1 102 
+2054 (bottom panel) and SDSS 0922+2928 (top panel). The uncertainties in the fluxes from 
the SDSS photometry are comparable to the size of the dots. In both cases, the w-band flux 
from the photometry lies below the extrapolation of the spectroscopy. The best-fit models 
to the photometry for an assumed logg = 8.0 are shown as blue lines (T e g = 8270 K for 
SDSS 0922+2928 and T cS = 8150 K for SDSS 1 102+2054; see Tables 2 and 3). Shown as red 
lines are the best-fit models to the absorption lines for an assumed log g = 8.0 (T cff = 10 000 K 
for SDSS 0922+2928 and T cff = 10 500 K for SDSS 1102+2054). As part of the fitting the 
absorption lines, the models were normalised to the observed SDSS spectra. Here, we show 
the best-fit models scaled to the r-band magnitude without that normalisation to illustrate the 
systematic uncertainty in the optical slope/effective temperature. 



16 




6200 6400 6600 7700 7800 7900 8200 




4000 5000 6000 7000 8000 9000 

Wavelength [A] 



Supporting Figure 2. Model fits to the observed absorption lines of SDSS 1 102+2054, illustrat- 
ing the effect of varying the surface gravity. The central panel shows the full wavelength range 
of the SDSS spectroscopy, with a model for (T e s = 10 500 K, \ogg = 8.0) and the correspond- 
ing abundances from Table 3 superimposed on the SDSS data. Note that the continuum slope is 
normalised to the observed data, and that only the absorption lines are fitted. Offset by two and 
four flux units are spectra with (T cff = 11 500 K, log g = 8.5) and (T cff = 12 500 K, log g = 9.0), 
respectively, keeping the abundances fixed to those given in Table 3. The top panels show 
close-ups of the normalised SDSS spectrum around several O I and C I lines, as well as of Ha. 
Line identifications are indicated by tickmarks. A model for (T cS = 10 500 K, log g = 8.0) 
is superimposed on the observed spectrum, models with (T cS = 11500K, log g = 8.5) and 
(T eff = 12 500 K, log g = 9.0) are offset by -0.4 and -0.6 flux units, respectively. 
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